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ABSTRACT 

The oxygen abundance gradients among nebular emission line regions in spiral galaxies have been 
used as important constraints for models of chemical evolution. We present the largest-ever full- 
wavelength optical spectroscopic sample of emission line nebulae in a spiral galaxy (M31). We have 
collected spectra of 253 HII regions and 407 planetary nebulae with the Hectospec multi-fiber spec- 
trograph of the MMT. We measure the line-of-sight extinction for 199 HII regions and 333 PNe; we 
derive oxygen abundance directly, based on the electron temperature, for 51 PNe; and we use strong 
line methods to estimate oxygen abundance for 192 HII regions and nitrogen abundance for 52 HII 
regions. The relatively shallow oxygen abundance gradient of the more extended HII regions in our 
sample is generally in agreement with the result of iZaritskv et a l. (1994), based on only 19 M31 HII 
regions, but varies with the strong-line diagnostic employed. Our large sample size demonstrates 
that there is significant intrinsic scatter around this abundance gradient, as much as <~ 3 times the 
systematic uncertainty in the strong line diagnostics. The intrinsic scatter is similar in the nitrogen 
abundances, although the gradient is significantly steeper. On small scales (dcprojcctcd distance < 0.5 
kpc), HII regions exhibit local variations in oxygen abundance that are larger than 0.3 dex in 33% 
of neighboring pairs. We do not identify a significant oxygen abundance gradient among PNe, but 
we do find a significant gradient in the [N II] ratio that varies systematically with surface brightness. 
Our results underscore the complex and inhomogeneous nature of the ISM of M31, and our dataset 
illustrates systematic effects relevant to future studies of the metallicity gradients in nearby spiral 
galaxies. 

Subject headings: Catalogs — Galaxies: individual (M 31) — Galaxies: evolution — Galaxies: abun- 
dances — HII regions: general — Planetary nebulae: general 



1. INTRODUCTION 

The galactocentric radial gradient of chemical abun- 
dance within spiral galaxies has become an important 
parameter in modeling th e chemical evolution of galaxies 
(|Henry fe Wort hev 1999]). The gradient is the manifesta- 
tion of a variety of physical processes acting from galaxy 
formation to the present, including gas infall, star for- 
mation history, stellar initial mass function, and radial 
migration. Historically, the observation of abundance 
gradients has motivated the development of an alytic 
models of chemical evolution (jLvnden-Belll I1975L and 
references therein) an d served as a valuable constraint 
for d et ailed modeling (iMaciel et al~l 12005b iMolla fc Di'azl 
2001 iScannapieco et al.l 120081: iMagrini et al.l 120091: 



Cariei fc Peimbertl 120100 . Moreover, the characteriza- 



tion of the abundance profiles of spiral galaxies has sig- 
nificant applications to other proble ms, including the in- 
terpr etation of stellar populations (iMassev fe Johnsonl 
Il998fl and the stud y of supernova explosion sites 
(|Levesaue et al.ll201ll ). 

The oxygen abundance of gas throughout the disk of 
star-forming galaxies can be measured via optical spec- 
troscopy. O is the most easily accessible metallicity proxy 

nsandcrs@cfa.harvard.edu 



due to its high relative abundance and its strong opti- 
cal emission lines from both its major ionization states 
(0 + and ++ ). A variety of diagnostics have been de- 
veloped to estimate O abundance from the f lux ratios 
of prominent optical emission lines (see e.g. |Stasiriska 
2002tlKewlev fe Ellisodl2008l: lLopez-Sanchez fe Esteban 
20101) . These include "direct" methods, whereby the 
electron temperature of the nebula is derived from 
measurement of the weak auroral line [O III] A4363 
(|Osterbrock fe Ferlandl [20061: [G arne ttl [1992T1 . and " sta- 
tistical" or "strong line" measures, where the abundance 
is inferred from ratios of only the brightest elemental 
lines. Different abundance diagnostics carry systematic 
discrepancies as large as 0.7 dex (a factor of 5), which 
must be carefully consi dered when interpreting results 
(|Kewlev fc Ellison! [20081) . 

Strong line methods have enabled extensive obser- 
vational studies of the abundance distribution of dis- 
tant spiral galaxies. While high quality spectra can 
be obtained for the nearby HII regions of the Milky 
Way, galactocentric distance can be determined more 
easily for extern al galaxies (iStanghellini et aU 120081 : 



IHenrv et al.l[20Tol) . fehields fc SearTd (119781) derived a 



negative (outward-decreasing) abundance gradient from 
just three HII regions in M101, and by this time it 
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was already expected from theory that negative radial 
abundan ce gradients would be characteristic of all spiral 
galaxies. iZaritskv et al.l (|1994l) compiled radial O abun- 
dance profiles for 39 spiral galaxies, finding gradients 
ranging from to —0.23 dex kpc -1 — only the abundance 
profile of NGC 2541 was inconsistent with a negative or 
flat gradient. 

Because PNe abundances should reflect the older ISM 
of their progenitors, differences between the HII re- 
gion and PNe abundance gradients can be used to in- 
fer time-vari ation in the radial abundance trend of a 
galaxy (e.g. IMagrini et~aT] 120091 for M33). A sample 
of high-luminosity PNe should reflect current oxygen 
abundances in a galaxy, by selecting massive progenitors 
which would have formed recently. The oxygen abun- 
dances in a sample containing only the brightest PNe 
in a galaxy would represent populations with ages from 
~3x 10 7 — 1 10 years earlier tha n the HII regions in the 
same galaxy ([Magrini et al.|[2009() . For surveys penetrat- 
ing more deeply, the proportion of PNe from less massive 
progenitors will grow quickly both because less massive 
stars are more common and because their resulting PNe 
are longer lasting. However, stars of mass < 1.5 M Q may 
not form PNe because their envelop e is ejected so s lowly 
as to disperse before being ionized (Stasihska 2002). 

While ~ 10 3 bright HII regions (Lh q > 5 x 
10 34 e rgs s -1 ) a re kno w n in M31 dBaade fe Arp| 



1961 IPellet et all Mfik IWalterbos fc Braunl 11992: 
Azimlu et al.l 12011 ). previous spectroscopic surveys to 



determine abundance have provided a bunda nce esti- 
mates for < 7% of them. iBlair et "all (|1982T ) found a 
significant radial abundance gradient in a sample of 11 
HII regions: O/H decreases by a factor of 4 from about 
4 — 23 kpc in galactocentric distance. This gradient is in 
ag reement with an earl i er stu dy of 8 HII regions in M31 
bv lDennefeld fc Kunth (119811) . In a n analysis of the data 
from both Blair et al. (Il982f) and iDennefeld fc Kunthl 
(|1981h . IZaritskv et al] (|1994[) quantified the oxygen 
abundance gradient as relatively shallow among their 
sample of galaxies: —0.020 ± 0.00 7 dex/kpcB- In com - 
parison, the largest gradient in the [Zaritsky et al.l (|1994D 
sample was -0.231 ± 022 de x/kpc (NGC 3344). Most 
recently, iGalarza et al.l (|1999t) surveyed 46 HII regions 
out to about 20 kpc in the Northeastern portion of 
M31, find ing a radia l gradi ent in R23 that is consis- 
tent with IBlair et "ail (|1982D . but only among objects 
wh ose morphology was classifi ed as " center- brightened" 
by IWalterbos fc Braunl (|1992D . Among HII regions of 
other morphologies, they report a flat abundance pro- 
file. Each of these previous surveys has relied on var- 
ious strong line abundance diagnostics. Additionally, 
abundance gradients have been derived f rom s urveys of 
M31's stellar population. iTrundle et"al~l (|2002fl found a 
flat radial metallicity profile among seven young B stars. 
iWorthev et al.l ()2005l ) estimated abundances for an older 
population, red giant stars, from color-magnitude dia- 
grams of 11 fields in M31. For radii < 25 kpc, they 
find a negative grad i ent si milar to the nebu l ar res ult of 
IDennefeld fc Kunthl (|l98l and IBlair et al.l (fl982l ). but 

1 Herea fter, we quote physical values for M31 based on a dis tance 
of 770 kpc HFreedman fc MadorellT990t) . IZaritskv efaLI lu"994T ) used 
an earlier distance measurement of 700 kpc, so we apply a 10% 
correction to their measured abundance gradient. 



they report a flattening at larger radii. 

M31 hosts as ma ny as 10 4 PNe, about tw ice as many 
as the Milky Way (iNolthenius fc Fordlll987D . In a kine- 
matic survey, IMerrett et al.l (|2006l) cataloged 2615 of 
M31's PNe. Abundances have been previously derived 
for less than 1% of these. IJacobv fc Fordl (|1986| ) deter- 
mined abundances for three PNe: two in the halo and 
one in the outer disk. The results of the largest previous 
spectroscopic survey of M31 P Ne, including 70 obj ects, 
have not yet been published (jKniazev et al.l I2005h . A 
survey of 30 nebulae in the bu lge of the galaxy was per- 
formed bv lRicher et al.l (|1999l ). but only a lower limit of 
abundance could be derived for 14 of the nebulae. Be- 
cause the survey only extended out to < 4 kpc and be- 
cause only the bulge population is sampled. iRicher et al.l 
(|1999D di d not investigate the r adial abundance gradient. 
Recently, IKwitter et al.l (|2012|) have derived abundances 
for 16 PNe in the outer disk of M31. 

By providing a more thorough characterization of the 
abundance profile of in M31, this paper seeks to enable 
an improved understanding of the chemical evolution of 
M31 and similar spiral galaxies. In Tables [HQ] we present 
the largest available spectroscopic catalogs of HII regions 
and PNe in M31. Our observations come from the Hec- 
tospec multifibcr spectrograph on the MMT, whose mul- 
tiplexing ability provides a large advantage over previous 
surveys of these objects. In Section[2]we describe the ob- 
servational parameters and analytical techniques used to 
produce the catalog. In Section[3]we compare our results 
to previous publications and discuss trends and implica- 
tions identified in the catalog. Our principal findings are 
that a significant negative abundance gradient is only 
demonstrated among the brightest and most diffuse HII 
regions. The radial profile is more flat among dimmer 
or more compact HII regions, and in general there is a 
large amount of scatter in the physical properties of the 
ISM of M31. We characterize this scatter in terms of the 
radial distribution of extinction and abundance in HII 
regions and PNe, and also in terms of the discrepancies 
among neighboring objects. We provide a summary of 
major results in Section 2J 

2. OBSERVATIONS 
2.1. Data collection 

Small, resolved objects and unresolved Ha features 
were selected as HII region candidates by inspecting 
the images of the L ocal Group Galaxies Survey (LGGS, 
iMassev et al.ll2007T ). Additionally, some objects observed 
as part of the M31 cluster survey of lCaldwell et al.l (|2009l ) 
had strong emission, and are included in the present 
st udy. Some objects fr om the planetary nebula catalog 
of IMerrett et al.l (|2006l) as well as strong and unresolved 
[O III] features from the LGGS images were observed as 
PN candidates. We have excluded from o ur sample any 
object s identified as emission line stars by IMassev et al.l 
(|2007t) . one object which was found to have broad emis- 
sion lines characteristic of supernova remnants, and ~ 20 
objects which showed broad emission features character- 
istic of WR stars. 

We note that many M31 HII region s have diameters 
greater than 50" (|Arp fc Bru cckcl 1971), far too large to 
be encompassed by a single Hectospec fiber, which sub- 
tends 1.5"on the sky (~ 5 pc at the distance of M31). 
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Moreover, to facilitate a future study of the kinematic 
properties of nebulae in M31, we attempted to avoid HII 
regions which would have a large internal velocity dis- 
persion. For both these reasons, the largest HII regions 
were therefore intentionally omitted from the sample, al- 
though some objects re main (mostly the s tar clusters) 
that were estimated by iBlair et aTT (|1982[ ) as at least 
48 pc in diameter. We assume that any inhomogeneities 
in spectral properties within each HII region are small 
and random and therefore that the region sampled by 
the Hectospec fiber is representative of the whole (see 
e.g. IMcCall et all 119851 IPellegrini etaHIMnh . 

Optical spectra were obtained with the Hectospec 
multi-fiber posit ioner and spectrogra ph on the 6.5m 
MMT telescope (jFabricant et all 120051) . The Hectospec 
270 gpm grating was used and provided spectral cover- 
age from 3650 — 9200A at a resolution of ~ 5A. Some 
spectra did not cover [O II]A3727, because of the design 
of the spectrograph (alternate fibers are shifted by 30A) , 
and the small blueshift of M31. The observations were 
made in the period from 2004 — 2011 as a component 
of a survey previously described in lCaldwell et all ([2009D 
and were reduced in the uniform manner outlined there. 
The frames were first debiased and flat fielded. Indi- 
vidual spectra were then extracted and wavelength cal- 
ibrated. Standard star spectra obtained intermittently 
were used for flux calibration and instrumental response. 
Sky subtraction is achieved with Hectospec by averaging 
spectra from "blank sky" fibers from the same exposures 
or by offsetting the telescope by a few arcseconds (see 
ICaldwell et aL|[2009D . Because local background subtrac- 
tion could potentially subtract object flux from extended 
HII regions, we compare the spectra reduced from re- 
peated observations of the same objects using different 
sky spectra, local and distant. For both stellar and dif- 
fuse objects, we find only small differences in the ratios of 
Ha/H/3 and also [N II] /Ha. The rms of the discrepancy 
in the log of the flux ratio is ~ 0.05 dex and the mean 
difference is ~ 0.01 dex. We therefore conclude that the 
sky subtraction is adequate. 

Each of 25 1 degree fields in M31 was exposed for be- 
tween 1800 and 4800 s. The spectra of objects that were 
observed multiple times (in overlapping fields) have been 
combined, effectively summing those integration times. 
Sample spectra are shown in Figure [TJ The locations of 
emission nebulae included in this study are shown in Fig- 
ure [2j Many of our HII regions fall in the "Ring of Fire," 
a circular feature visible in HI density m aps extending 
from about 8 to 15 kpc in the disk of M31 (jSofue k. Katol 

[Mi- 
Velocities were measured using the SAO xcsao soft- 
ware and emission line templates (one typical of HII re- 
gions and another of PNe). Repeat measurements of 
114 objects gave an rms of a single measurement of 2.1 
k m s -1 . We also comp ared our velocities with the work 
of IHallidav et al | (I2006D . who also u sed a multifibcr spec- 
trograph. and IMerrett et"aTI (|2006l ) who used the Plane- 
tary Nebula Spectrog raph (which imaged in [O III] ) . The 
IHallidav et all (|2006f ) comparison revealed a mean offset 
of +7.1 km s " 1 wit h an rms of 5.4 km s _1 , while the 
iMerrett et all (|2006f ) comparison resulted in a mean off- 
set of +2.9 km s -1 with an rms of 15.6 km s _1 . The 
IMerrett et all (|2006f ) rms was expected to be larger be- 
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Fig. 1. — Three spectra from Hectospec. From the bottom, a 
high-metallicity HII region (object HII166), a low- metallicity HII 
region (object HII153), and a PNe (object PN130). The two HII 
regions are separated by only 1.93' on the sky, corresponding to 
a separation of ~ 0.4 kpc at the distance of M31. For clarity, the 
relative flux scales were set such that the strongest emission line 
in each spectrum (Ha for HII regions and [O III]A5007 for the PN) 
are the same height and the wavelength zero point was offset by 
100 and 200A for the two upper spectra. [O II]A3727 falls off of 
the original spectrum for the upper and lower spectra. 

cause of their use of the single emission line. These com- 
parisons led us to estimate our mean error to be 3 km s _1 . 
The final step in the reduction was the transformation of 
the spectra to zero velocity using the observed velocities. 

2.2. Classification 

These candidates were formally classified according 
to two-dimensional emission line ratio tests of excita - 
tion mechanism (BPT diagrams, iBaldwin et "all 1 1981 1 ) . 
Specifically, we apply equation 1 of (|Kniazev et al.ll2008|) . 
which distinguishes HII regions from PNe based on 
their locations in a diagram of [O III]/H/3 (03) ver- 
sus [N II] /Ha (N2) (Figure [3]). These tests have the 
advantage of relying on ratios of strong emission lines 
that are near to each other in wavelength and therefore 
not sensitive to reddening corrections or instrumental re- 
sponse. The measurement of line fluxes is described in 
Section [231 We ado pt a slight amendmen t to the divid- 
ing line proposed by (jKniazev et alJl2~008l ) , illustrated by 
the dashed line in Figure [3l This amendment is favored 
because it classifies as PNe a number of objects that we 
have reason to believe are not HII regions. These ob- 
jects have stellar morphologies (see below) and many of 
them are in the halo, based on their position along the 
minor axis (\Y\ > 4 kpc), as is illustrated in Figure It 
is not surprising that there may be small offsets in the 
appropriate line ratio diagnostics for different samples 
given slight differences in flux-measurement methodol- 
ogy. The amended classification we adopt to distinguish 
HII regions from PNe is: 



03 > (0.61/(N2- 0.47)) + 1.0 



(1) 



In total, we identify 407 PNe and 253 HII regions 
among our spectroscopic sample. Of these objects, 
392 PNe and 200 HII region spectra were observed 
with the atmospher ic dispersion compensator (ADC, 
iFabricant etaTI [20081 ) . For some spectra, the ADC mal- 
functioned, and for those we do not present quantities 
derived from emission line flux ratios with large wave- 
length separations. 81 of the PNe do not appear to be 
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Fig. 2. — Locations of M31 HII regions and PNe observed with Hectospec, overlaid on the M31 mosaic from the Digitized Sky Survey 
(DSS). The HII region sy mbol s are colored to oxygen abundance, using the N06 N2 diagnostic, and the PNe are colored according to their 
N2 flux ratio (see Section 12,71 1. Halo PNe are denoted with squares and objects whose abundance/flux ratio could not be measured are 
marked in black. 

are such that it is not possible to derive their abundances, 
we remove them from the sample and do not discuss them 
further. An additional 25 objects could not be classified 
because the S/N in their [O III] and H/3 emission lines 
were too weak. We do not discuss these unclassified ob- 
jects further. 

In addition to this spectroscopic classification, we have 
made a morphological classification for each object based 
on its appearance in the LGGS Ha images. Objects 
whose Ha emission appears unresolved at the resolu- 
tion of these images (~ I") arc classified as "stellar" and 
those that appear extended and nebulous are classified 
as "diffuse." We note that this is similar t o the classifica- 
tion scheme of IWalterbos fc Braunl (|1992f ) referred to by 
IGalarza et al.l (|1999|) . where their type "C" or "center- 
brightened" corresponds to our "stellar" and their type 
"D" or "diffuse" corresponds to the same. We find that 
7 spectroscopically-classified PNe are diffuse. Because 
most of these objects are near the dividing line in the 
BPT diagram, it is difficult to determine whether these 
are PNe embedded in HII regions or simply poor spec- 
troscopic classifications. However, these only amount to 
~ 1% of our full sample of PNe. 

About half of our spectroscopically-classified HII re- 
gions are stellar. From Figure it is clear that the stellar 
HII regions are more likely to have high 03 or N2 ratios 
relative to diffuse HII regions and are therefore nearer 
to the dividing line in the BPT diagram. It is possible 
that some of these objects are PNe with an unusually 
low value of [O IIIl/H/? . Amo ng the sample of Galac- 
tic PNe in iHenrv et al.l (|2010[ ). there are some objects 
(Ml-11 and Ml-12) that fall on a similar region of the 
BPT diagram. Acc ording to the C atalogue of Galactic 
Planetary Nebulae (jKoh outek 2 0011) . thes e obje cts were 
originally classified as PNe by |Minkowsk]| (|1946t) , whose 
classifications were essentially morphological. It is also 
possible that these are extended HII regions whose Ha 
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Fig. 3. — Exc itation mechanism diagnostic diagram 

IIBaldwin et al.l 119811 ) distinguishing PNe from HII regions 
for M31 emission line nebulae observed in this study. Emission 
line flux ratio error bars are typically s maller than the points . The 
solid black classification divider is from [Kniazcv et al. (2003); the 
dashed line is our amended divider, on which our spectroscopic 
classification is based. The data points are color coded according 
to Ha morphology (stellar or diffuse); objects identified as being 
in the halo based on their position along the minor axis (\Y\ > 4 
kpc) are colored black. Objects shown on the edge of the figure 
denote flux ratio limits for lines which were not detected. 



projected onto the disk of M31 (according to the distance 
along the minor axis, \Y\ > 4 kpc, see Section l2~3j) . and 
we therefore assume they are associated with the halo 
population. 23 of the PNe are hosted by the dwarf satel- 
lite galaxy NGC 205 rather than in M31, according to 
their position and v elocity in the kinematic survey of 
iMerrett et al.l (|2006D . Because the S/N of these spectra 
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emission is simply dominated by the emission in the im- 
mediate vicinity of one bright star, and therefore appear 
compact on the LGSS images. Therefore the diffuse sub- 
set of HII regions in our sample may represent a cleaner 
sample, less likely to have PN contamination. 

2.3. Galactocentric distance 

The de projected gal actocentric distance was calculated 
followin g iHaudl (11981. We assume a distance to M31 of 
770 kpc dFreedman fc Madorelll990l ) and an inclination of 
i = 12.5° (jSimien et al.lll978fT . The adopted coordinates 
of the galactic cent er and posit ion angle of the major 
axis, precessed from IHaudl (|1981D . are: 

a = 00°42'44".52 (J2000) 

Sq = +41°16'08".69 (J2000) 

</>o = 37°42'54" 

Hereafter, we use "galactocentric distance" to refer to 
this deprojected distance. Any PNe with a distance pro- 
jected along the semiminor axis (Y) greater than 4 kpc, 
we associate with the halo population and do not con- 
sider in our sample when fitting for radial trends in the 
disk. Additionally, there are 3 objects that we consider 
halo PNe given their projected distance | V | > 5 kpc, de- 
spite the fact that they fall in the HII region regime of 
the BPT diagram. 

2.4. Emission line fluxes 

The line fluxes were measured using the line profile- 
fitting capabilities provided by the IRAF package 
f itprof s. We fit Gaussian profiles to a wavelength range 
20 A in width centered on the rest frame wavelength of 
each line. We fit a linear continuum to 20 A regions 
of the spectra off the wings of each line. For groups 
of nearby lines likely to be blended, we fit simultane- 
ous Gaussians. The line flux is estimated as the inte- 
gral of the fitted profile with the continuum subtracted. 
We estimate the uncertainty in the line fluxes using the 
Monte Carlo methodology implemented by f itprof s. In 
these Monte Carlo simulations, random noise is repeat- 
edly added to the spectrum according to a simple linear 
noise model dependent on two parameters: the gain of 
the CCD and a noise floor. We estimate the gain and 
noise floor for each individual line profile by compar- 
ing the spectrum of the background regions to the cor- 
responding variance spectrum. We have compared this 
estimate of the line flux uncertainty to the discrepancy 
between line flux ratios measured in repeated observa- 
tions of the same objects and find them to be compara- 
ble; for the Ha/H/3 ratio, the Monte Carlo uncertainty 
estimate and the discrepancy among the repeated obser- 
vations are each ~ 3% in the mean. Based on the Monte 
Carlo estimate of the uncertainty in the line flux, we re- 
quire that S/N> 3 and otherwise report a non-detection. 
Additionally, we record the equivalent width as measured 
by f itprof s, if the sky-subtracted continuum level is 
more than 2a greater than zero (where the continuum 
level and its standard deviation, a, are measured from 
the background regions of each line). 

The [O II] AA3726 - 3729 doublet is not resolved in 
these spectra, and the sum effectively measured will 



henceforth be referred to as [O II] A3727. Moreover, be- 
cause A3727 is on the very edge of our spectral range, it 
has a large flux-calibration uncertainty associated with 
it, which is propagated through to our estimate of the 
flux uncertainty. As mentioned above, for some spectra, 
A3727 is outside of the observed spectral range. 

We note that we have made no correction for underly- 
ing stellar absorption, but we estimate that this will not 
significantly affect the measurement of the Balmer line 
fluxe s. Using models from Starburst99 (jLeitherer et al.l 
Il999f ). we subtracted off underlying continua for popu- 
lations with ages ranging from 4 to 20 Myr for a rep- 
resentative low- and high-mctallicity HII region (objects 
HII166 and HII153), see Figure [1]). We then measured 
H7 and H/3. The ratios of those two lines changes at 
worst by 1% from the uncorrected values (the worst case 
results from assuming the youngest age for the underly- 
ing population). As a further test of the small effect that 
the underlying continuum has on our measurements, we 
also measured Ha equivalent widths and compared those 
with A v (derived from the Ha/H/3 emission line ratio), 
and found no significant correlation. This indicates in 
general that Balmer absorption doesn't significantly af- 
fect the measured Balmer ratios for the objects in our 
sample. 

We present the measured line fluxes and their uncer- 
tainties for each object in Table [2j For all parameters 
derived from these line fluxes, we propagate the uncer- 
tainty in the line flux via Monte Carlo simulations. In 
these simulations, we sample from a Gaussian line flux 
probability distribution with a mean and standard devi- 
ation as reported in Tabled] We then report the median 
and standard deviation of the resulting distribution as 
our best estimate of the derived parameters and its un- 
certainty. 

We characterize the surface brightness of each emission 
line region based on its observed Ha emission line flux, 
even if the objects are unresolved. Ha is used for this 
purpose because it is strong and easily detected in nearly 
every spectrum. However, the 1.5" Hcctospec fibers gen- 
erally do not cover the entire area of the diffuse emission 
nebulae. Moreover, as the observation nights were not all 
photometric, there could be photometric uncertainties of 
a factor of a few. We attempted to place all the spectra 
on the same photometric scale by using the multiple ob- 
servations, and find that this was successful to within a 
factor of 1.3. 

We therefore divide our sample into three surface 
brightness classes (1: "dim," 2: "normal", 3: "bright"; 
see Figure 0| based on the extinction-corrected Ha line 
flux, rather than assert a precise flux measurement. We 
establish separate brightness classes for the HII regions 
(based on the diffuse subset only) and PNe. For spectra 
where the H/3 line is not detected, we assume the me- 
dian Ay for that class (HII region or PNe) for the pur- 
pose of calculating the extinction corrected Ha line flux. 
The diffuse HII regions in our sample are at ~ 2.4 times 
higher surface brightness than the PNe, in the median (or 
~ 4.4 times brighter after extinction correction). More- 
over, the HII regions encompass a larger range in surface 
brightness, extending to nearly two orders of magnitude 
brighter in Ha. 

We have estimated the apparent brightness limits cor- 
responding to these surface brightness classes by match- 
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Fig. 4. — Top: The surface brightness distribution of diffuse HII 
regions (blue) and disk PNe (green), as measured spectroscopi- 
cally from the extinction-corrected Ho emission line flux (arbi- 
trary units). The solid and dashed vertical lines denote the surface 
brightness limits which split the sample into three equally sized 
bins of HII regions and PNe, respectively. Bottom: Comparison 
of spectroscopically-mcasured Het fluxes (not extinc tion corrected) 
and the photometric catalog of[Xzimlu et al. (20lJ) for HII regions 
only, with color coding by the catalog diameter. Symbol shape de- 
notes our morphological types. Points on the edge of the figure do 
not have matches within 5". The dashed line is the best fit among 
small (D < 20 pc) objects. 

ing the HII region s in our sample to the catalog of 
I Azimlu et aT1(|2011l ). using the nearest objects within 5". 
They use an automated code to segment diffuse emission 
in continuum-subtracted Haimages of M31 for the iden- 
tification of HII regions and PNe and to measure the Ha 
flux and diameter. Figure U shows the relation between 
our spectrosc opic Ha flux measur ements and the pho- 
tometry from lAzimlu et all (| 2 1 If ) . The scatter in the 
relation is correlated wi th the HII region d iameter, be- 
cause HII regions which lAzimlu et al.l (|2011l) segment as 
larger objects are have a correspondingly smaller frac- 
tion of their flux fall on our Hectospec fibers. Among 
small objects (D < 16 pc), the scatter is ~ 0.3 dex. 
This relation implies that our HII region flux bins have 



edges at 
[2.6,5.41 



[7.0, 31.1] x 10 
x 10~ 15 



ergs cm s and, for PNe, 



ergs cm s . 
2.5. Extinction 



A reddening correction was applied to restore the 
Balmcr recombination decrement of each spectrum to 
its theoretical value. We assume Ha/H/3=2.85, which 
corresponds to T=10,000K and n e — 10 4 cm -3 for Case 
B recombination. However, the Balmer ratios are not 
very sensitive to any of these parameters: the assump- 



tion of Case A only alters the ratio by < 0.05 and 
extreme temperatures and densities also have little ef- 
fect: Ha/H/3 = 3.04 for T = 5,000 K and n e = 



10 2 cm" 3 , Ha/H/3 = 2.73 for T 



20,000 K and 



n e = 10 6 cm' 3 (IQsterbrock fc Ferlaridl l2Q pj) . The ex- 
tinction curve given in lCardelli et al.l (J1989D was applied, 
with a value of 3.1 adopted for Ry. In cases where we 
derive a negative value of the extinction in the visual 
band (Ay < 0), we instead assume Ay = 0. 

The value of the extinction derived for each object is 
presented in Tables [3] and |U 

2.6. Direct abundance estimation 

We derive "direct" oxygen abundance estimates for 
PNe by estimating the electron temperature of the gas' 
dominant excitation zone, which can be done only if a 
temperature-sensitive line is detected. We use [O III] 
A4363 exclusively for this purpose; although auroral lines 
from other ions (e.g., [N II] A5755) are detected in a few 
spectra, their numbers are not sufficient for a statistical 
sample and their S/N is typically much lower than [O III] 
A4363. This prescription is also applied to our HII re- 
gions, however direct abundances can be derived for only 
4 HII regions due to the weakness of the [O III] A4363 
line. 

The electron temperatures are estimated using 
IRAF's five-leve l neb ular modeling package nebular 
(|Shaw fe Dufourl [19941) . The nebular task temden is 
first applied to iteratively estimate the ++ temperature 
(T e (0 ++ )) and density (n e ) of the nebula from the [O III] 
and [S II] line ratios, respectively When the [S II] lines 
are note available, we assume a reas onable range of dens i- 
ties, n e = 15±2 x 10 3 cm -3 (see e.g. lKwitter et al.ll2012D . 
If the measured line ratios correspond to unphysical con- 
ditions (outside the range for which temden is calibrated, 
500 < T e (0++) < 10 5 K and 1 < n e < 10 8 cm" 3 ), we 
do not calculate the direct abundance. The + temper- 
ature i s then es t imate d using the linear empirical rela- 
tion of lGarnettl (|1992t) . The 0+ and 0++ abundances 
are then estimated using the density, ionic temperatures, 
and [O II] and [O III] line ratios follow ing the ionizatio n 
correction factor (ICF) prescription of lShi et"aLl (|2006|) . 
The total oxygen abundance is taken to be the sum of 
these two ionic abundances. Th is methodology is sim ilar 
to that applied in, for example, IBresolin et al.l (|2010D . 

A variety of studies have shown that stellar evo- 
lution of PNe progenitors should not modify oxygen 
abundance at a level significant for th e identification of 
abun dance gradients from the nebulae (jRicher fc McCalll 
120071 and references therein). Furthermore, it has 
been shown that bright PNe (within ~ 2 mag of the 
brightest in the galaxy), presumably from more mas- 
sive stars with shorter lifetimes, have approximately 
the same oxygen abundances as the surrounding ISM 
and HII regions, to with in the observational uncertainty 
((Richer fc McCall|[2007l) . It has been demonstrated that 
O, as well as Ar and Ne, abundance gradients in the 
Milky Way as measur ed with HII regions are reflecte d 
in observations of PNe (jPottasch fc Berna rd- S alas 20061 ). 
However, oxygen may be dredged up in some low- 
mctallicity cases (log(0/H)+12 < 8) and thereby enrich 
the PNe relative to its progenitor by Alog(Q/H)< 0.3 
dex (jRicher fc McCalll 120071: iHernandez- Martinez et al.l 
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120091: iMagrini fe Goncalvesl I2009T) . In high-mass PNc, 
oxygen may be depleted by Alog(0/H)< 0.1 dex dur- 
ing the lifetime of the progen itor via the ON-cyclc 
(jHernandez-Martinez et aT1 l2009'). Moreover, it has been 
observed that the N abundanc e of PNe will often ex- 
ceed that of the local ISM tiRicher fe McCalll [20071 
iHernandez-Martinez et aTl I2009T ) — justified theoretically 
by nitrogen production via the CN or ON cycle d uring 
the first and second dredge ups. iHenrv et all (poTot) have 
established empirically that the oxygen abundance gradi- 
ent of Milky Way PNe does not depend on the Peimbert 
type of the PNe or the vertical distance from the plane 
of the Galaxy; we therefore do not consider these param- 
eters here. 

2.7. Strong line diagnostics 

HII region abundances were calculated from the 
extinction-corrected line flux ratios according to a num- 
ber of independently calibrated abundance diagnos- 
tics from the literature. These diagnostics each de- 
pend on a different combination of emission line ra- 
tios. Each calibration has its own characteristic scatter 
and systematic offset as compared to the other methods 
(jKewiev fe Ellisonll2008T ). It is unclear that any particu- 
lar method is superior, and some methods may be more 
applicable than others for certain comparisons. In partic- 
ular, diagnostics tied to the direct abundance scale may 
agree better with ste llar oxygen abundan ces within the 
same galaxy (see e.g. IBresolin et al][20 09;). We therefore 
employ multiple methods and keep their differences in 
mind when discussing results. We report the abundance 
derived from each method in Tables [3] (HII regions) and 
H] (PNe). 

First, we apply the R 23 oxygen abundance calibration 
of lZaritskv et all (|1994| ). an average of three earlier meth- 
ods, hereafter referred to as "Z94." Z94 is only calibrated 
for the higher-metallicity upper branch of the well-known 
i?23-abundance degeneracy. The majority of M31 HII 
regions may be expected to fall on this upper branch, 
given that all the M3 1 HII regions in the compilation of 
iZaritskv et al.l (|1994| ) did. If the measured line ratios cor- 
respond to an abundance outside of the range for which 
Z94 is calibrated (generously, 8.4 <log(0/H)+12 < 9.6), 
we do not record the measurement. 

Second, we apply the [N Hi /JO II] oxygen abundance 
calibration of lKewlev fe Dopital (120021) . her eafter referred 
to as "KD02." iKewlev fe Dopital (|2002T ) synthesize a 
variety of modern photoionization models and observa- 
tional calibrations to produce recommendations for pro- 
ducing an abundance estimate given different permuta- 
tions of available emission lines. We i mplement the pre- 
script ion outlined in the appendix of lKewlev fc Ellison! 
(j2008| ). as follows. We use the [N 11} /[O II] ratio 
to break the degeneracy between the upper and lower 
branches of R23 ■ For the upper branch, we emplo y the 
[N II] /[O II] calibration of lKewlev fc Dopital ([200l . For 
the rare l ower branc h case s, we average the R23 diag- 
nostic s of iMcGaughl (|1991| ) and iKobulnickv fc Kewlevl 
(|2004[ ). If the measured line ratios correspond to an 
abundance outside of the range for which KD02 is cali- 
brated (8.2 <log(0/H)+12 < 9.6), we do not record the 
measurement. 

Third, we apply the empirical [N II] /Ha ("N2") and 
[O III]/[N II] ("03N2") oxygen abundance calibrations 



of iNagao et~aT1 (|2006l ). hereafter referred to as "N06." 
We prefer the N06 diag nostic to the simi l ar "P P04" N2 
and 03N2 diagnostics of iPettini fc Pagell (|2004[ ) because 
N06 is well-calibrated in the high-metal licity regime of 
M31 using data from the SDSS galaxies (jTremonti et al.1 
120041 ) . If the measured line ratios correspond to an abun- 
dance outside of the range for which N06 is calibrated 
(7.0 <log(0/H)+12 < 9.5), we do not record the mea- 
surement. Because this diagnostic relies on the N II lines 
to measure the O abundance , scatter is introduced by 
varia tions in the N/O ratio (jPerez-Montero fc Continl 
120091 ). We note also that the 03N2 diagnostic is not re- 
liable when 03N2 > 100, due to line saturation, but this 
only occurs in a metallicity regim e lower than that sam- 
pled here (log(Q/ H)+12 < 7.5; IPettini fc Page! l200l 
INagao et al]l2006l ). 

Fourth, we apply the excitation parameter 
("P method") ox ygen abundance calibration of 
iPilvugin fc Thuanl (|2005D . hereafter referred to as 
"PT05." This is an updated version of the calibra- 
tion first defined in IPilvugin! (|2001bf ) (P01). P is 
calculated from the ratio of [O III] to ([O II] + [0 III]). 
PT05 additionally relies on the R23 line ratio, so the 
[N II] /[O II] ratio is used to break the R23 degeneracy. 
If the measured line ratios correspond to an abundance 
outside of the range for which PT05 is calibrated 
(6.8 <log(0/H)+12 < 9.1), we do not record the 
measurement. 

Fifth, we ap p ly the nitrogen abundance calibration of 
IPilvugin et"all (j2010f ) hereafter referred to as 'PVT," 
relying on the combination of P and [O III], [N II], 
and [S II] ("ONS") line ratios. The PVT diagnostic is 
calibrated separately for each of three different [N II] 
regimes. If the measured line ratios correspond to an 
abundance outside of the range for which PVT is cali- 
brated (7.3 <log(N/H) + 12 < 8.9), we do not record the 
measurement. 

Wc have not attempted to factor in the systematic er- 
ror in the abundance diagnostics, although they are typ- 
ically much larger (<~ 0.1 dex) than our reported errors, 
which are derived by pr opagation of the l i ne flu x un- 
certainties. For example, IKewlev fc~ Ellison (2008) esti- 
mates the rms scatter between relative metallicities mea- 
sured with the Z94 diagnostic, as compared to other pop- 
ular diagnostics, is 0.07 dex based on a sample of 30,000 
SDSS galaxies. They find that the scatter in the other 
diagnostics are similar, the largest mean rms belonging 
to P01 (related to PT05) at 0.11 dex. Because the only 
references for the accuracy of each abundance estimation 
technique are estimates from other diagnostics, which are 
not necessarily independent, any quantification of uncer- 
tainty must be interpreted with caution. 

We have compared the oxygen abundance measure- 
ments made for the same HII region in different diagnos- 
tics. Among the strong line methods, there is very good 
agreement between the Z94 and KD02 methods (stan- 
dard deviation of 0.07 dex). There is fairly good agree- 
ment between Z94 and N06 N2 (median offset of -0.10 
dex and standard deviation of 0.21 dex) and between the 
N06 N2 and 03N2 diagnostics (negligible median offset, 
standard deviation of 0.12 dex). The PT05 diagnostic 
does not agree well with the other strong line methods, 
having a median offset as large as -0.50 dex (Z94) and a 
standard deviation as large as 0.30 dex (N06 N2). 
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In Figure [5] we show the cumulative distribution func- 
tions (CDFs) of oxygen abundance for the HII regions 
and PNe in our sample as derived by the different diag- 
nostics. The total range in PNe abundances is about 
7.6 < k>g(0/H) + 12 < 8.8 (a factor of 16). The 
range of HII region abundances varies widely by di- 
agnostic. For example, PT05 abundances range from 
8.0 < log(0/H) + 12 < 8.5 while N06 N2 abundances 
range from 7.9 < log(0/H) + 12 < 9.5. 

The difference in both the shape and median value of 
the CDFs of different diagnostics are due to two factors: 
the systematic discrepancy between the calibrations, and 
the selection effects imposed by the requirement for cer- 
tain emission lines to be detected in order to apply 
each diagnostic. One additional selection effect is the 
range over which the diagnostics arc calibrated; for exam- 
ple, the Z94 abundance scale is only calibrated down to 
log(0/H)+12 = 8. 4, as described above. If t he diagnostic 
transformations of iKewlev fe Ellison! (|2008l ) are applied 
(along with the trivial transformation between PP04 N2 
and N06 N2), the CDFs for the N06 N2, KD02, and Z94 
diagnostics agree reasonably well (Kolmogorov-Smirnov 
p- value > 0.01), despite the selection effects. The dis- 
crepancy is largest between the PT05 diagnostic and the 
other strong line methods. The median abundance for 
the N06 N2 diagnostic is log(0/H)+12 = 8.89, while for 
PT05 it is 8.34 dex. While this offset is approximately 
equal to the offset determined for the SPSS galaxies with 
these diagnostics (jKewlev fc Ellisonl 120081 and also con- 
sidering the conversion from N06 N2 to PP04 N 2 and 
PT05 to P01). However. IKewlev fc Ellisonl (|200l) deter- 
mined that no reliable transformation can be established 
between PT05 and the other diagnostics because the rela- 
tion is highly non-linear and has high scatter (> 0.1 dex). 
From Figure [5] it is clear that the direct abundances may 
only be derived for the lowest-metallicity objects where 
the auroral line is accessible. 

The combination of selection effects and calibration 
discrepancies should be kept in mind as we investigate 
the oxygen abundance profile of M31. Results should 
only be compared if they are quoted in the same diag- 
nostic, due to the systematic deviations and scatter that 
exist between different diagnostics. Moreover, if a diag- 
nostic does not probe the full range of abundances in the 
population, then the measured abundance profile in the 
diagnostic should not be used to infer physical proper- 
ties of the galaxy. In particular, we suggest that care be 
used in interpreting the abundance profile derived from 
the direct diagnostic, because it imposes severe selection 
effects (the availability of the auroral line) that selec- 
tively exclude high-metallicity objects from samples, and 
the PT05 diagnostic, because it is produces abundances 
discrepant with other strong line diagnostics. However, 
measurements from these diagnostics are still useful. The 
direct abundance is the only diagnostic based on a mea- 
surement of the true electron temperature of the nebulae, 
and it provides a us eful comparison to Galactic work (e.g. 
IShaver et al.l Il983h . Moreover, the direct method must 
be applied to PNe, for which strong line methods are not 
available. The PT05 diagnostic is the most modern cali- 
bration of the P method, which typically reflects the oxy- 
gen abundance as measured by the dire ct method more 
closely than other strong line methods (lPilvuginll2001al 
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Fig. 5. — The cumulative distribution functions (CDFs) of oxygen 
abundance for M31 HII regions and PNe (both disk and halo) as 
derived by different diagnostics. The strong line methods are shown 
as solid lines, while the direct abundances are dashed. The number 
(N) of objects for each diagnostic is noted in the legend and is 
limited by the available line flux ratios from each spectrum and 
the abundance range over which the diagnostic is calibrated. PNe 
include halo objects. A color version of this figure is available in 
the electronic journal. 

as well as our results above). 

3. DISCUSSION 

Our analyses are primarily concerned with looking for 
radial trends in the interstellar medium properties of 
M31, with the goal of identifying any information that 
describes the chemical evolution history of the galaxy. 
For this purpose, we focus on objects in the disk of the 
galaxy. While all HII regions studied in this survey are 
attributed to the disk, a large population (N = 81) 
of PNe appear in projection to be outside of the disk 
(\Y\ > 4 kpc); this is the halo population discussed in 
Section l2~3l We exclude these disk PNe from our analy- 
sis, except where explicitly described. 

We look for radial trends in the optical extinction, oxy- 
gen abundance, and nitrogen abundance of HII regions 
and PNe in M31. Our initial analyses, fitting linear 
trends and looking for correlations, are summarized in 
Table [5] This table summarizes the significance of radial 
trends in two different ways: 

Bootstrap: We fit a line by ordinary least squares to 
the radial distribution of the parameter, and then 
repeat many times with resampling. We simultane- 
ously resamplc from the set of all objects with mea- 
surements of that parameter (with replacement) 
and also from the probability distribution function 
of the derived parameter. In this way, we can es- 
timate the slope and intercept (extrapolated value 
at the center of M31) of the radial trend in a way 
that is not sensitive to outliers or objects with poor- 
quality spectra. 

Spearman: We report the Spearman rank correlation 
statistic p and its p value. A p value much less 
than zero indicates a strongly negative gradient. A 
p value much greater than zero would indicate that 
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any apparent correlation with radius could be due 
merely to chance. 

We discuss these gradient analyses in the following sec- 
tion. 

3.1. Extinction 

The estimated values of extinction in the visual band 
(Ay) as derived from the Balmer decrement are reported 
in Tables [3 ] and [H and vary from to nearly 5 mag. 
iBlair et alT(|1982ft s ampled HII re gions with extinction 
up to ~ 2.3 mag and lGalarza et aLl ()1999h up to ~ 4 mag. 

Figure |6] illustrates the Ay distribution versus galac- 
tocentric distance for 199 HII regions and 333 PNe. The 
extinction is patchy. Because reddening imposes a selec- 
tion effect on our sample, the maximum extinction we 
observe at a given radius is just a lower limit. Nonethe- 
less, the maximum extinction varies radially: as high as 
~ 5 mag from ~ 10—15 kpc, in the Ring of M31, to 
< 1.5 mag beyond 20 kpc, in the outskirts of the disk. 
However, objects fill the figure down to ~ mag at all 
radii. The maximum extinction of PNe follow a similar 
radial trend, but typically have smaller extinction than 
HII regions ( Ay,PNe < 3 mag) . It is expected that HII re- 
gions will have higher extinction than PNe, because they 
arc near the large dust clouds associated with star forma- 
tion and because they are primarily found in spiral arms 
rather than evenly throughout the disk. PNe from the 
halo population are also likely to be found above the disk 
of M31 where extinction should be negligible. Moreover, 
HII regions are typically an order of magnitude brighter 
than PNe a nd therefore m a y be observed t hrough greater 
extinction UPanagial 1 1978ft . iKumarl (|1979f ) asserted that 
a plateau at < 1 mag in the extinction of HII regions 
begins at 12 kpc, while our data demonstrate that large 
(> 2 mag) values of extinction are common out to nearly 
twice this distance. 

PNe in the halo population have consistently small 
extinction values. The median (and 16th, 84 per- 
centile) value for extinction among halo PNe is A^haio = 
0.08^08 mag for AT halo = 75, while for disk PNe it is 
much larger: Ay disk = 0.47lo 4g mag for N disk = 333. 

It is reasonable to expect the extinction to trace 
abundance, because dust grains that cause redden - 
ing are composed of heavy elements (jShieldsl I1990D . 
Some studies h ave reported s uch parallel gradients in 
spiral galaxies (ISarazi iJUzi IViallefond fe Gossl Il98l 
Ivan der Hulst et al.lll988l . for M33. M101, and M51, re- 
spectivel y), while others have m e asured flat extinction 
profiles (IViallefond fc Gossl [19861 Kaufman et "all H987L 
for M33 and M81). In Figure [7] we show the oxygen 
abundance of HII regions and PNe in M31 against the 
extinction (Ay). This figure illustrates that there is not 
a clear correlation between extinction and oxygen abun- 
dance among the objects in our survey. 

Regardless of any trends in the extinction with radius 
or abundance, it is clear that the extinction in M31 is 
patchy. For example, neighboring HII regions (< 0.5 
kpc in deprojected distance, or « 2.2 ' separation on the 
sky) differ in extinction by as much as 2.9 mag. Among 
the 98 such neighboring pairs in our sample with extinc- 
tion measurements, 33% have a discrepancy in Ay of 
more than 1.1 mag (Figure [TT]) . Some of this deviation 
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Fig. 6. — Radial extinction profile for HII regions and disk PNe 
in M31. The procedure for estimating Ay is described in the text. 
Typical error bars are smaller than the size of the points, but 90th 
percentile error bars are shown to illustrate the uncertainty for the 
least-constrained PNe (green) and HII regions (blue). 
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Fig. 7. — The oxygen abundance of HII regions (N06 N2) and 
disk PNe (direct) in M31 vs the extinction (Ay) as measured from 
the Balmer decrement. Statistical error bars (derived as described 
in the text) are shown for each dimension, but are typically too 
small to be visible. 

is attributable to the large inclination of M31 introduc- 
ing somewhat disparate column densities into the line of 
sight for apparently adjacent objects that are in front of 
and behind the disk. 

3.2. Radial oxygen abundance gradient 

The measurement of any abundance gradient among 
the M31 HII regions depends strongly on the choice 
of abundance diagnostic. For example, if we em- 
ploy the N06 N2 diagnostic, we measure a gradient 
(-0.0195 ± 0.0055 dcx kpc" 1 ) that is negative at the 
~ 4<r-level and consistent wit h the canonica l valu e 
(-0.020 ± 0.007 dcx kpc- 1 ) of IZaritskv eTail (fl99l . 
If we instead use the Z94 diagnostic, we find a gradi- 
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ent that is much less steep and only different from zero 
at the - la level (-0.0208 ± 0.0048 dex kpc -1 ). The 
N06 03N2 (-0.0130 ± 0.0068 dex kpc" 1 ) and KD02 
(—0.0096 ± 0.0049 dex kpc -1 ) diagnostics yield simi- 
lar results. The p- value of the Spearman test sug- 
gests that a real correlation exists in all four cases, to 
varying degrees (p N2 = 3 x 10 -04 , pz94 = 2 x 10 -06 , 
P03N2 = 0.06, Pkdo2 = 1 x 10 -04 ). If the PT05 diagnos- 
tic is used, however, we do not find a significant gradient 
(-0.0054 ± 0.0064 dex kpc" 1 ). The Spearman test re- 
flects the lesser significance of the radial trend in this 
diagnostics (pptos = 0.60). Moreover, different results 
can be achieved if the sample is divided by morphologi- 
cal type or surface brightness (Section 13. 4p . We do not 
detect the temperature-sensitive auroral lines for enough 
HII regions to investigate the abundance gradient in the 
direct diagnostic. 

An illustrative radial oxygen abundance profile of M31 
HII regions is shown in Figure [51 using the N06 N2 di- 
agnostic. This diagnostic is highlighted because it relies 
on only the brightest emission lines and is insensitive 
to flux calibration and reddening correction. It there- 
fore produces reliable abundance estimates for a very 
large number of HII regions (N — 192). Also shown is 
the abundance gradient as fit by the bootstrap method: 
(8.98 ±0.08) dex + (-0.0195 ± 0.0055) dex kpc" 1 . This 
represents a relatively shal ow gradient am o ng th e nearby 
spiral galaxies studied by iZaritskv et al.l (|1994|) , falling 
in the 11 — 26 percentile range based on the la er- 
rorbars quoted above. Employing the isophotal radius 
po = 16 kpc they define for M31, the size normalized 
gradient falls in the 16 — 32 percentile range. 

The radial oxygen abundance profile of M31 PNe is 
shown in Figure [HI using the direct abundance diagnos- 
tic. Also shown in the figure is the abundance gradient 
as fit by the bootstrap method, (—0.0056 ± 0.0076) dex 
kpc -1 , which is consistent with zero. The Spearman p- 
value is fairly large (0.45), emphasizing that there is no 
significant correlation between the PNe abundances and 
galactocentric radius. However, various systematic ef- 
fects influence the interpretation of this correlation, as 
we will discuss in Section [3.51 

3.3. Intrinsic scatter 

It is clear from Figures [8J and [9] that, regardless of what 
the true slope of the abundance gradients may be, there 
is significant intrinsic scatter about the trend. 

In Figure QJJ] we characterize this scatter by calculat- 
ing the standard deviation of the abundance of all HII 
regions in different radial bins for each abundance diag- 
nostic. By dividing the sample into radial bins, as op- 
posed to calculating the standard deviation of the entire 
sample, we partially remove the variance that would be 
introduced by an abundance gradient. While the KD02, 
PT05, and Z94 diagnostics seem to produce the least 
scatter (~ 0.1 dex), this could be in part due to se- 
lection effects; the emission lines necessary for calculat- 
ing the R23 ratio are not accessible in fainter HII re- 
gions. If we consider the diagnostics with the largest 
sample size, N06 N2, the scatter in abundance rises from 
~ 0.2—0.3 dex from the inner to outer regions of the disk. 
This is significantly larger than the scatter inherent to 
the abundance diagnostics themselves (e.g. ~ 0.07 dex 



for Z94 and PP04 N2. lKewlev fc Ellison! [20081 ). 

The conclusion that the intrinsic scatter in the abun- 
dance gradient is larger than the observational uncer- 
tainty reflects some studies in the literature. If the best- 
fit gradient is subtr acted from the M31 H II region abun- 
dances measured in IZaritskv et al.l ([1994D . then the stan- 
dard deviation among the abundances is ~ 0.16 dex — 
similar to what we measu re with the Z94 diagnostic. 
iRosolowskv k, Simonl ([2008D found an intrinsic scatter of 
0.11 dex among their 61 HII regions in M33, and as- 
serted that this is lar ger than t he pr ecision of the mea- 
surements. However. IBresolinl (|2011l ) have argued that 
certain systematic ef fects have artificially increas ed the 
scatter measured by IRosolowskv fc Simonl (|2008f ) , par- 
ticularly the inclusion of high-excitation HII regions and 
low S /N spectra. Because only 4 of our spec t ra me et the 
strict S/N threshold suggested by IBresolinl (|2011|) (e.g. 
S/N ([O III] A4363) > 5), our dataset is not sufficient 
to address the intrinsic scatter in a subset of the data as 
they recommend. However, any high-e xcitation o bjects 
in our sample of the type discussed by IBresolinl ([201 If ) 
would instead be classified as PNe (Section I2.2p . and 
therefore would not contaminate the HII region statistics. 
Moreover, the uncertainties in the strong line abundance 
measurements are negligible compared to the measures 
intrinsic scatter. For example, for the N06 N2 diagnostic 
our measurement uncertainties as propagated from the 
emission line flux uncertainties have a median of 0.03 dex 
and 90th percentile value 0.10 dex, much smaller than 
the > 0.2 dex intrinsic scatter we measure in the radial 
abundance profile. Future spectroscopic studies of HII 
regions in M31 should acquire spectra of sufficient S/N 
in [O III]) A4363 to measure the intrinsic scatter in the 
direct diagnostics. 

In Figure [UJ we investigate local fluctuations in the 
ISM of M31. We do so by considering the discrepancy in 
extinction and abundance measurements among HII re- 
gions and disk PNe separated by less than w 2.2 ' on the 
sky, corresponding to < 0.5 kpc in deprojectcd distance 
in M31. As we have previously discussed in Section |3~T1 
local fluctuations in extinction are often quite large — 
A(0/H)> 0.3 dex and AA V > 1.1 mag for one third of 
HII regions and A(0/H)> 0.2 dex for one third of disk 
PNe. 

As illustrated in Figure [2] the oxygen abundance of 
the ISM of M31 is inhomogencous. Neighboring HII re- 
gions (< 0.5 kpc in deprojected distance) differ in oxygen 
abundance (N06 N2) by as much as 0.6 dex, an order of 
magnitude. Among the 132 such neighboring pairs in 
our sample with N06 N2 abundance measurements, 33% 
have a discrepancy in log(0/H) of more than 0.3 dex 
(Figure ITT]). These discrepancies could be partially ex- 
plained by measurement uncertainty, however 0.4 dex is 
~ 5 x the scatter expected from the systematic uncertain- 
ties in the dia gnostic (~ 0.07 dex for the similar PP04 
N2 diagnostic. iKewlev fc Ellism] |20M) . This scatter is 
similar to the maximum discrepancy among the 8 PNe 
and 5 HII regions (~ 0.3 and ~ 0.2 dex, respectively) 
in the immediate solar neighborhood (< 2 kpc) observed 
bv lRodriguez fc Delgado-Ingladal (|201lt ) . Only 20 neigh- 
boring PNe pairs have direct abundance estimates, so we 
do not consider their distribution of discrepancies. 

We present spectra for an example of two neighbor- 
ing HII regions with discrepant abundances in Figure [TJ 
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Fig. 8. — The radial oxygen abundance profile of HII reg ions in M31 using the N06 N2 diagnostic. Each point is color-coded by the surface 
brightness of the Ha emission line, as discussed in Section 12.41 The thick line is the best-fit linear abundance gradient, as derived using the 
bootstrap method described in the text. The lower and upper dashed lines are constructed from the 33rd and 66th percentile slope and 
offset parameters calculated in the bootst rap simulation, respec tively, shown to illustrate the possible variation in the fitted gradient. The 
dotted line is the canonical gradient from [Zaritsky ct al. (1994). Statistical error bars (derived as described in the text) are shown for the 
oxygen abundance, but arc frequently too small to be visible. 



These are two diffuse HII regions (objects HII153 and 
HII166) that are separated by only 1.93 ' on the sky, cor- 
responding to a separation of ~ 0.4 kpc at the distance 
of M31. When we calculate their galactocentric radii, 
the difference is 0.84 kpc, and their velocities only differ 
by 24 km s^ 1 . Despite being so nearby, object HII153 is 
low metallicity (log(O/H) N06 N 2 + 12 = 8.36 ± 0.04) and 
object HII166 is high metallicity (log(O/H)N06 N2 + 12 = 
8.94 ±0.02). 

While it would be interesting to search for non-linearity 
in the abundance profile (such as breaks near the well- 
known star-forming ring of M31), the high-level of intrin- 
sic abundance scatter present throughout the disc would 
make it difficult to evaluate different models. 

For PNe, we estimate the intrinsic scatter in the O 
abundance of M31 disk PNe as > 0.26 dex. Because 
we find no evidence of a significant radial trend in oxy- 



gen abundance for PNe, we simply calculate this number 
as the standard deviation of the 51 disk PNe with di- 
rectly measured abundances. We consider this to be a 
lower limit because high-mctallicity PNe are systemat- 
ically excluded from our direct abundance sample due 
to the weakness of the auroral line. The median oxygen 
abundance is log(0/H)+12 = 8.46 dex. For the 17 PNe 
in the halo of M31 for which we can measure direct abun- 
dances, we find a median value and standard deviation of 
log(0/H)+12 = 8.50 ±0.22 dex. The median abundance 
in the halo is therefore lower than that in the disk, but 
this distinction is small given the intrinsic scatter in the 
a bundances of each p opulation. 

iHenrv et all (|2010l ) derive the oxygen abundance gra- 
dient in the Galaxy from observations of 124 PNe 
with high-quality spectra and well-determined distances. 
They report a best-fit gradient of —0.058 ± 0.006 dex 
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Fig. 9. — The radial oxygen abundance profile of PNe in M31 
using the direct diagnostic. The sample o f PN e is divided into 3 
luminosity classes as described in Section 12.41 A radial gradient 
is fit to the sample of each flux class, and the dashed red line 
illustrates the gradient fit to the full sample. 
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Fig. 10. — The scatter in the radial oxygen abundance profile 
of the diffuse HII regions in M31. The horizontal dashed line at 
0.07 dex represents the typical systematic error in the diagnostic 
(the mean rms error in the Z94 diagnostic from[Kcwlcy fc Fllisonl 
120081) . The scatter is the standard deviation of the abundance of 
all HII regions in each of three different radial bins. The bins are 
chosen such that each bin has an equal number of objects. The 
bin edges are therefore different for each diagnostic, because the 
sample changes when different line ratios are required. Various 
abundance diagnostics are employed, with sample sizes as noted in 
the legend. 

kpc -1 . They assert that the scatter around this best-fit 
gradient is ~ 40% larger than the uncertainties they as- 
cribe to the abundance estimates. Similarly, wc find an 
intrinsic scatter in PNe abundances that is larger than 
the measurement error, the scatter in the M31 disk PNe 
abundances we report above is > 2x larger than the 
median uncertainty in our disk PNe direct abundance 
estimates (a ~ 0.10 dex, as derived by propagation of 
the line flux uncertainties) . An accounting of systematic 



errors could inflate the asserted measurement error, but 
they could not account for the intrinsic scatter unless 
t hey alter the measured abun danc es by a factor of > 2. 

iRosolowsky fc Simonl (|2008t ) and IMagrini et all (|20T0f ) 
invoke inefficient azimuthal mixing to explain local fluc- 
tuations in th e ISM metallicity of M33 (but see also 
lBresolinll2011| ). In this scenario, mixing performed by ve- 
locity shear due to differential rotation occurs on a longer 
timescale (~ 10 s yrs) than enrichment by star formation 
in the spiral arms of the galaxy. Such a scenario could 
also apply to the inhomogeneities we observe in the ISM 
ofM31. 

3.4. Dependence on HII region properties 

Due to its size and proximity, M31 provides a unique 
laboratory for studying the ISM of a spiral galaxy; for 
many extragalactic studies, only the brightest HII re- 
gions in the galaxy are accessible to spectroscopy. While 
the large number of relatively-dim HII regions included 
in our survey (Section |2.4|) allows us to probe the ISM 
properties of M31 more thoroughly than ever before, it 
also has the potential to introduce discrepancies with 
past work. Here we investigate whether the measured 
abundance profile varies systematically with the bright- 
ness or compactness of the HII regions in the sample. 

In Figure [12] we investigate the fitted abundance gra- 
dient parameters (slope and characteristic abundance at 
12 kpc) for the diffuse HII regions as a function of Ha 
emission line flux density by dividing our sample into 
surface brightness bins (see Section 12 4|) . Although the 
bins are comprised of equal numbers of HII regions, be- 
cause the S/N of spectral lines depends strongly on sur- 
face brightness, there are typically fewer abundance mea- 
surements available in the lower surface brightness bins. 
Wc adopt a minimum of 5 abundance measurements for 
performing abundance gradient analysis, which is the 
mi nimum number of HII regions sampled for any galaxy 
bv lZaritskv et all (|1994f ). There are not sufficient PT05 
abundance measurements to perform this analysis with 
that diagnostic. We find that the fitted slope and char- 
acteristic abundance parameters are essentially indepen- 
dent of surface brightness, varying by amounts consistent 
with the error bars on the fitted parameters. 

This resul t cont rasts with that reported by 
Ma grini et all (|2010f ) . who found that the abun- 
dance gradient in M33 was more than twice as steep for 
bright "giant" HII regions than for the less-luminous 
objects in their sample. Given the distance of M33 
(840 kpc), their threshold for bright HII regions 
{Fft a ~ 1.2 x 10 34 ergs s _1 ) is similar to our threshold 
for normal/bright HII regions. They attribute the 
surface-brightness dependence of the gr adient to self- 
enrich ment in giant HII regions. However, Magrini et al. 
(|2010f ) derive metallicity gradients based on direc t abun- 
dance measurements. Because IMagrini et al.l (|2010f ) 
exclude ~ 1/3 of their spectroscopic sample due to 
insufficient S /N in the emission lines necessary to derive 
direct abundances, one possible systematic explanation 
for the stronger gradient they measure among brighter 
HII regions is that the [O III] A4363 line is weaker in 
high-metallicity HII regions. Therefore high-mctallicity 
HII regions may be systematically excluded from 
their lower-luminosity (not giant) HII region sample. 
Similarly, the strong-line methods wc apply in this 
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Fig. 11. — The distribution of discrepancies between the extinction value and oxygen abundance (N06 N2 for HII regions, direct for PNe) 
of neighboring HII regions and PNe. Neighbors are defined as any pair of objects that are separated by < 2.2' on the sky, corresponding 
to < 0.5 kpc in deprojected distance in M31. Halo PNe are excluded. There are 98 such neighboring HII region pairs with extinction 
measurements and 132 such pairs with oxygen abundance measurements in our sample; 502 and 20 PNe pairs, respectively. 



study carry certain selection effects (see Section 12. 7|) . 
but we do not find that the gradient parameters vary 
significantly with brightness. 

We can also investigate the abundance gradient as a 
function of HII region morphology (as described in Sec- 
tion [22]). As described in Section VTM the morpholog- 
ical and surface brightness classifications are quite dif- 
ferent, and t he stellar HII region s may have PN con- 
tamination. iGalarza et al.l (|1999f ) only found a signif- 
icant radial abundance trend in the R23 line ratio for 
"center-brightened" (stellar) HII regions, and no gradi- 
ent for HII regions of more extended morphology (dif- 
fuse). Generally, the abundance gradient slopes for 
diffuse HII regions reported in Table El are in good 
agreem ent with the canonical value froni lZaritskv et al.l 
(|1994f ). while the slopes for stellar HII regions are more 
shallow or consistent with zero. For example, in the 
Z94 diagnostic, the Spearman test indicates a signif- 



icant gradient among diffuse HII regions (p = 6 x 
1(T 05 ), but not stellar (p = 0.06). In detail, Tabled 
demonstrate that for most abundance diagnostics (e.g. 
N06 N2), the radial gradient slopes derived from ob- 
jects of either morphology agree to within their error- 
bars (Alog(O/H) + 12 N06 N2 [diffuscstollar] = [-0.0195 ± 

0.0070 , -0.0224 ± 0.0082] dex kpc" 1 ). IGalarza et all 
(|1999f ) do not report a trend in [N II] /Ha for objects of 
any classification, whereas we report significant and simi- 
lar abundance gradients for this line ratio in both diffuse 
(AlogQN II]/Ha)diftuso =-0.0177 ± 0.0049 dex kpc" 1 ) 
and stellar (AlogQN II]/Ha) ste iiar =-0.0080±0.0056 dex 
kpc" 1 ) HII regions, equivalent to the gradient measured 
in N06 N2 oxygen abundance. We do find a systematic 
offset in the N06N2 diagnostic, such that the best fit gra- 
dient intercept (metallicity at the galactic center) for the 
[diffuse,stellar] HII regions is log(0/H) + 12 =[9.09±0.09, 
9.24 ± 0.11] in the N06N2 diagnostic. This can be ex- 
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Fig. 12. — The fitted parameters of the oxygen abundance 
(log(0/H)+12) gradient for diffuse HI1 regions in M31 in bins of 
Ho; emission line flux density. The bins were chosen such that they 
each have an equal sample size, as described in Section 12.41 The 
dotted vertical lines denote the width of the bin. The points are 
given small x-offscts for clarity. The y-errorbars are the bootstrap 
parameter uncertainties as described in the text. Different abun- 
dance diagnostics are used as noted in the legend. 

plained by a systematic difference in the hardness of the 
ionizing radiation in HII regions of these two morpholog- 
ical classes, which causes strong line methods which do 
not account for the ionization parameter to overestimate 
the metallicity of compact HII regions. 

However, there are particularly strong disagreements 
between stellar and diffuse HII regions among the KD02 
and PT05 diagnostics (Table [5]). In detail, these discrep- 
ancies are driven by the presence of anomalously high- 
metallicity (log(0/H) +12 > 9.0 on the KD02 scale) stel- 
lar HII regions. Because these HII regions are also dim, 
with surface brightness consistent with PNe in our sam- 
ple, they are likely to be PN contaminants. 

3.5. Time-variation in the abundance gradient 

The HII regions in our sample are in general more en- 
riched than the PNe, as is demonstrated by Figure 
Similarly, the difference in the a bundances of HII re - 
gions and PNe in M33 reported by|Magrini et al.l ([20100 
was ~ 0.1 dex and interpreted in the context of the 
time-varying composition of the ISM. The median oxy- 
gen abundance and standard deviation for M31 disk PNe 
is log(0/H)+12 = 8.46 ± 0.26 dex. Among HII regions 
studied with the N06 N2 diagnostic, the typical abun- 
dances are much larger (log(0/H)+12 = 8.89 ± 0.24). 
This would suggest a similar discrepancy between HII 
regions and PNe as in M33; however, the median abun- 
dance we measure for PNe may be depressed because of 
the selection effect on the [O III] A4363 line required to 
estimate direct abundances. Moreover, using a diagnos- 
tic that typically correlates between with direct abun- 
dance measurements, PT05, we find a smaller median 
metallicity for HII regions (log(0/H)+12 = 8.34 ±0.12); 
again, selection effects should act to exclude highcr- 
metallicity objects. This illustrates the complicating role 
of systematic effects in comparing abundance measure- 
ments for a statistical sample of extragalactic HII regions 



and PNe. 

In Figure [131 we bin the PNe by surface brightness to 
investigate the potential time-dependence of the abun- 
dance gradient. As discussed in SectionQ] the luminosity 
of PNe are related to the masses of their progenitor stars 
and therefore to their ages. As in Section [331 wc adopt 
a minimum of 5 measurements for performing gradient 
analysis in each bin, and we are therefore not able to fit 
an abundance gradient for the least bright PNe. 

If the radial abundance gradient has strengthened 
over time, we would expect to find a gradient that 
is more negative with increasing PN surface bright- 
ness. In fact, the trend in the direct abundance is 
never strongly inconsistent with zero for either surface 
brightness bin. Among the brightest PNe, the slope 
is -0.0023 ± 0.0097 dex kpc" 1 . We do find that the 
median metallicity increases with brightness class, with 
[8.26,8.42,8.49] dex for [Dim, Normal, Bright] PNe. How- 
ever, because the auroral lines are only detectable in 
lower-metallicity PNe, selection effects could be elimi- 
nating high-metallicity PNe from our sample which could 
reveal a significant abundance gradient. Moreover, this 
selection effect would act more strongly to remove high- 
metallicity objects among the dimmer PNe, which would 
mimic the signature of increasing metallicity with PN 
surface brightness. 

Because direct abundance estimates can only be made 
for 51 of our 326 disk PNe spectra, it is worthwhile to 
look for gradients in strong line ratios such as R23 and 
[N II] /Ha. While these PN line ratios are not directly 
abundance-sensitive as they are in HII regions, significant 
gradients (if present) could indicate interesting trends in 
other physical properties. For [N II] /Ha, the gradient is 
positive at the ~ 2a level for the lowest-surface bright- 
ness bin (0.0178 ± 0.0108 dex kpc -1 ), and increases sig- 
nificantly with surface brightness such that it is positive 
among the brightest PNe (-0.0122 ± 0.0059 dex kpc" 1 ). 
The existence of this N2 gradient, and its correlation with 
surface brightness, could indicate a time-varying gradi- 
ent in excitation, chemical composition, or both. For 
R23, the slope is never significantly different from zero in 
any brightness bin. 

It would be difficult to distinguish time evolution in 
the abundance profile of M31 by comparing its PNe to 
HII regions, due to the intrinsic scatter in the popula- 
tions and the uncertainties in the determination of the 
abundance gradients. Results would be particularly in- 
fluenced by the choice of strong line diagnostic for HII 
region abundances, and by the cuts made on morphol- 
ogy/surface brightness (Section [ 



3.6. HII region nitrogen abundance gradient 

While the oxygen abundance gradient is the most 
observationally accessible (Section [l}, the radial nitro- 
gen gradient is of particular interest because it may be 
the steepest of any observable el ement (e.g. in M33 , 
iMagrini et alll2010t i. The models of lMagrini etaTl (|2010f ) 
show that the N gradient should be steeper than the O 
gradient due to the different timescales for production; 
N is produced primarily in low and intermediate-mass 
stars, while high-mass stars more efficiently produce O. 

The radial nitrogen abundance profile of M31 HII 
re gions is shown i n Fig ure [T4l using the diagnostics 
of iPilvugin et al.l (|2010l ). Also shown is the highly- 
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Fig. 13. — The fitted parameters of the oxygen abundance 
(log(0/H)+12) gradient and line ratios for PNe in M31 in bins of 
Ha emiss ion line flux density. The figure is constructed similarity 
to Figure fl2l 

significant (~ 4a) radial gradient for the diffuse objects 
as fit by the bootstrap method, (7.83 ± 0.08) dex + 
(-0.0300 ± 0.0058) dex kpc -1 . This slope is approxi- 
mately as steep or steeper than all the oxygen abundance 
gradients reported for any diagnostic and morphological 
selection in Table [5] However, it is only one third as 
steep as the nitrogen gradient in M33 (—0.08 ± 0.03 dex 
kpc- 1 . IMagrini et alj l2010h . If we include the stellar 
HII regions (N = 10), we find a gradient that is simi- 
lar (-0.0303 ± 0.0049 dex kpc" 1 ). 

As for the oxygen abundances, we find a large intrin- 
sic scatter about this gradient. Subtracting the fitted 
trend among all the HII regions, we find an rms scat- 
ter of 0.11 dex. This is significantly larger (<~ 2 times) 
the systematic uncertainty attribut able to the strong line 
diagnostic of (jPilvugin et al.l 12010ft . who report an rms 
scatter of 0.05 dex (smaller than for the equivalent oxy- 
gen abundance diagnostic). 

3.7. Comparison to previous observations of M31 

A variety of authors have previously derived abun- 
dance gradients for M31 from other surveys of HII re- 
gions, as well as O and B s tars. 

As we h av e disc uss ed. IDennefeld fc Kunthl (|1981l ). 
Blair et all (11982ft iZaritskv etHI (|1994j ). and 
Galarza et al.l (|1999ft have previously derived the 
abundance gradient of M31 from surveys of HII regions. 
From those works, we have the cano nical result of 
-0.02 ± 0.007 dex kpc -1 , produced by IZaritskv etlll 
(fl99l fr om the combined sampl e of 19 HII regions from 
IDennefeld fc Kunthl (fl98l and IBlair et all (fl982ft. In 
their survey of 46 HII regions in M31. iGalarza et al.l 
(|1999ft found agreement with this gradient, but only 
among center-brightened HII regions. Among HII 
regions of other morphologies, they found no significant 
tr ends. 

iTrundle et all (l200l derived oxygen abundances for 
7 stars of type O and B from about 5 to 30 kpc in 
the disk of M31. Additiona lly, they provided a re- 
analysis of the HII regions of IBlair et all (fl982ft using 



Fig. 14. — The radial nitro gen abundance profile of HII regions 
in M31 using the diagnostic of Pilyugin et al. (2010). The figure is 
constructed similarly to Figure [8] 

modern diagnostics, producing gradients that range from 
-0.013 (P-method) to -0.027 (M91) dex kpc" 1 . Their 
least squares fit to the stars yielded an abundance gra- 
dient of — 0.017 ± 0.02 dex kpc -1 , in good agreement 
with the canonical result. However, they noted that 
if they omitted the possible multiple system OB 8-76, 
then they instead measured a slope consiste n t with zero 
(-0.006 ± 0.02 dex kpc" 1 ). iWorthev et all ((2005ft pro- 
duced color-magnitude diagrams from HST images to de- 
rive the stellar abundance gradient in M31, finding (with 
some uncertainty in their zero point) good agreement be- 
tween the upper bound of the metallicity o f the stellar 
popul ation and the HII regions surve yed by IBlair et al.1 
i| 1982ft and IDennefeld fc Kunthl (|1981ft at a given radius. 

The results of previous surveys of HII regions and high- 
mass stars suggest that any discrepancy between the M31 
abundance gradients derived from different sources are 
dependent on systematic effects related to sample se - 
lection and diagnostic. Similarly, lUrbaneia et al.l (|2005f ) 
have found that the nebular abundance gradient can be 
smaller than, similar to, or larger than the stellar abun- 
dance gradient in the late-type spiral NGC 300 depend- 
ing on the choice of diagnostic. 

4. CONCLUSIONS 

We have presented optical spectroscopy of an unprece- 
dented sample of HII regions and PNe in a massive spi- 
ral galaxy, M31. In total, we reported line flux mea- 
surements for 253 HII regions and 407 PNe. We have 
derived the extinction, nitrogen abundance, and oxygen 
abundance for subsets of these objects using a variety of 
methods, as described in Section [21 From the analysis of 
these observations, we emphasize the following conclu- 
sions: 

1. For HII regions, we find an oxygen abundance 
gra dient generally consi stent with that found 
by IZaritskv et"afl (|1994ft . Using the N06 N2 
diagnostic, we find a gradient of (—0.0195 ± 
0.0055) dex kpc -1 among 100HII regions. We find 
a significantly steeper gradient in nitrogen abun- 
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dance, (-0.0303±0.0049) dex kpc" 1 among 52 HII 
regions. These represent relatively shallow gradi- 
ents as compared to other nearby spiral galaxies. 

2. For PNe, we detect no significant oxygen abun- 
dance gradient among 51 objects for which the [O 
III] A4363 line is detected and the direct method 
can be applied. However, using the line ratio [N 
II] /Ha, which is measurable for most PNe (N = 
277), we find significant gradients which vary sys- 
tematically with PN brightness, as illustrated in 
Figure [TU 

3. The ISM of M31 is highly inhomogencous. Both 
the visual extinction (Ay) and oxygen abundance 
vary significantly among even very nearby HII re- 
gions (2.9 mag and 0.6 dex for some HII regions 
separated by < 2 kpc, see Figures l2l and [TTj) . More- 
over, the intrinsic scatter observed about the HII 
region oxygen abundance gradient (~ 0.1 — 0.3 dex, 
see Figure I10[) is larger than the uncertainty in- 
herent to the strong line diagnostics (~ 0.1 dex). 
Similarly, the scatter among PNe in our sample is 
> 0.26 dex. 

4. The abundance gradient derived for HII regions in 
M31 is dependent upon the strong line metallicity 
diagnostic employed, and can be affected systemat- 
ically by sample characteristics such as HII region 



morphology and surface brightness. In particular, 
for observations to a given depth, some diagnostics 
can only be applied to low or high-metallicity HII 
regions, unless they have sufficient surface bright- 
ness. Among more compact (not extended, i.e. 
"stellar") nebulae that are spectroscopically clas- 
sified as HII regions, we find evidence of PN con- 
tamination that can lead to erroneous strong line 
abundance measurements. Sample characteristics 
of this type should be taken into careful considera- 
tion to mitigate systematic effects in future surveys 
of the ISM of nearby spiral galaxies. 
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TABLE 1 

Basic data for M31 HII regions and PNe 



ID RA DEC R (kpc) a Morph. Type b SB C ADC d Velocity (km s" 1 ) M06° RBC f AMB g 



HII regions 



HII001 


0:37:24.12 


+40:17:56.2 


23.0 


s 


2 


n 


-483.6 


HII002 


0:37:29.91 


+40:15:37.2 


21.9 


s 


1 


n 


-517.5 


HII003 


0:37:47.35 


+39:51:30.8 


23.7 


s 


1 


n 


-492.6 


HII004 


0:37:59.17 


+40:15:37.2 


19.2 


s 


1 


n 


-476.2 M2372 


HII005 


0:38:22.51 


+40:10:52.8 


18.4 


s 


2 


n 


-564.5 


HII006 


0:38:39.79 


+40:34:48.0 


18.2 


d 


1 


y 


-469.3 


HII007 


0:38:41.29 


+39:47:40.0 


28.2 


s 


1 


y 


-526.8 


HII008 


0:39:03.75 


+39:53:28.9 


26.9 


s 


1 


y 


-479.0 


HII009 


0:39:07.69 


+40:40:05.1 


16.2 


s 


2 


y 


-486.1 


HII010 


0:39:13.09 


+40:41:13.9 


15.9 


s 


1 


y 


-488.5 



HII9 



Hill 
Hill 



Planetary nebulae 



PN001 


0:38:44.18 


+40:17:58.8 


16.6 


1 


y 


-465.4 


M2364 


PN7 


PN002 


0:39:02.52 


+40:22:50.5 


15.2 


3 


y 


-529.4 


M2964 


PN9 


PN003 


0:39:05.33 


+40:41:45.9 


17.1 


2 


y 


-387.9 


M3198 




PN004 


0:39:06.59 


+40:14:59.7 


17.1 


3 


y 


-560.0 


M2357 




PN005 


0:39:14.80 


+40:24:28.5 


14.6 


1 


y 


-338.0 


M2304 


PN1 


PN006 


0:39:15.00 


+40:26:34.5 


14.3 


3 


y 


-538.8 


M2308 


PN1 


PN007 


0:39:15.80 


+40:12:38.4 


18.3 


2 


y 


-488.8 


M2972 




PN008 


0:39:16.29 


+40:22:12.8 


15.0 


2 


y 


-551.0 


M2960 


PN1 


PN009 


0:39:18.51 


+40:09:18.9 


19.9 


1 


y 


-379.4 


M2404 




PN010 


0:39:26.12 


+40:44:25.9 


15.2 


1 


y 


-407.7 


M1967 





Planetary nebulae (halo population) 



PNhOOl 


0:35:50.74 


+42:21:04.4 


104.8 


2 


y 


-296.8 


M2543 


PNh002 


0:36:27.15 


+42:06:21.8 


89.9 


1 


y 


-15.2 


M2549 


PNh003 


0:37:09.28 


+42:38:18.9 


103.6 


1 


y 


-245.0 


M2542 


PNh004 


0:37:28.34 


+42:10:57.6 


83.3 


1 


y 


-231.0 


M2544 


PNh005 


0:37:54.37 


+42:14:49.3 


81.7 


1 


y 


-334.2 


M7 


PNh006 


0:38:48.46 


+41:39:37.3 


51.2 


2 


y 


-195.4 


M352 


PNh007 


0:38:56.63 


+39:47:14.1 


29.7 


2 


y 


-417.0 




PNh008 


0:39:01.08 


+41:51:11.2 


56.5 


3 


y 


-216.0 


M174 


PNh009 


0:39:21.16 


+42:25:16.3 


75.2 


1 


y 


-331.8 




PNhOlO 


0:39:31.59 


+42:11:56.7 


65.0 


1 


y 


-343.5 


M5 



Unclassified 



X001 


0:39:06.10 


+40:37:22.8 


15.7 


s 


1 


y 


-512.5 


Hill 


X002 


0:39:07.30 


+40:36:25.5 


15.4 


s 


1 


y 


-513.9 


Hill 


X003 


0:39:08.89 


+40:24:11.5 


14.8 


s 


1 


y 


-534.1 


Hill 


X004 


0:39:09.40 


+40:29:16.0 


14.4 


s 


1 


y 


-542.2 


Hill 


X005 


0:39:12.90 


+40:50:59.6 


19.7 


s 


1 


y 


-408.2 


Hill 


X006 


0:39:16.09 


+40:43:55.5 


16.4 


s 


1 


y 


-530.5 


Hill 


X007 


0:39:16.90 


+40:20:55.5 


15.4 


s 


1 


n 


-560.4 


Hill 


X008 


0:39:24.30 


+40:48:07.5 


16.8 


s 


1 


y 


-469.8 


HII2 


X009 


0:39:34.80 


+40:49:59.5 


16.1 


s 


1 


y 


-446.3 


HII2 


X010 


0:40:10.29 


+40:45:19.0 


10.4 


s 


1 


y 


-533.0 


HII5 



Note. — The tabic is divided into sections based on the spectroscopic classification of the object as either an HII region 
(HII), planetary nebula (PN), PN in the halo population (PNh), or unclassified (X). This classification is based on line ratio di- 
agnostics, as described in §2.2\ Table ^ is publi shed in its ent irety in the electronic edition of the Astrophysical Journal and at 
https://www.cfa.harvard.edu/-nsariders/papers/M31/summary.html A portion is shown here for guidance regarding its form and content. 
a The galactoccntric radius of the object, calculated as described in ^2.31 

The morphological classification of the object based on LGGS Ha imaging, as described in ^2~2] 

The surface brightness class of the object as defined in fl2.4| bascd on the Ha flux. 

The status of the atmo 
indicates that it was not. 

c The ID number of the object from the PNe catalog of lMerrett et all J2006D . 
^ The name of the object from version 3.5 of the Revised Bologna Catalog, 
e The ID number of the object from the HII region catalog of |Azimlu et a! 



Galleti et all J2007I) . 
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TABLE 2 

Measured line fluxes for M31 HII regions and PNe, relative to H/3 



ID 



[O II] 
A 3727 



[O III] 

A 4363 



[O III] 
A 4959 



[O III] 

A 5007 



[N II] 
A 6548 



Ha 
A 6562 



[N II] 

A 6584 



[S II] 
A 6717 



[S II] 
A 6731 



HII regions 



HII001 


50 ± 20 








48 ±3 


442 ±5 


150 ±3 


57 ±4 


45 ±5 


HII002 










54 ±6 


255 ±4 


151 ± 5 






HII003 










26 ±6 


264 ±6 


76 ±7 


22 ±4 


20 ±4 


HII004 






40 ± 10 


100 ± 10 




252 ±8 


70 ±8 






HII005 










87 ±9 


409 ±9 


270 ± 10 






HII006 


430 ± 40 






110 ± 20 




340 ± 20 








HII007 






299 ±6 


852 ±9 




213 ±6 








HII008 




15 ±4 


286 ±5 


830 ± 10 




221 ± 5 








HII009 








50 ± 10 




606 ± 10 








HII010 


270 ± 10 








27 ±6 


306 ±6 


82 ±6 


43 ±6 


33 ±6 


Planetary nebulae 


PN001 






383 ±5 


1115 ± 10 




247 ±9 








PN002 






462 ±5 


1360 ± 10 


29 ±3 


411 ± 5 


75 ± 3 






PN003 






149 ±4 


442 ±4 




312 ±5 


23 ±6 






PN004 


34 ±3 


13 ± 1 


440 ±4 


1300 ± 10 


16 ±2 


292 ±3 


47 ±2 






PN005 






470 ± 10 


1410 ± 10 




320 ± 10 








PN006 


70 ± 20 




312 ±5 


930 ± 10 


50 ±6 


443 ±6 


146 ±6 






PN007 


110 ± 20 




191 ±3 


562 ±6 


39 ±8 


322 ±8 


113 ±8 






PN008 


80 ± 10 




381 ±5 


1160 ±10 


27 ±7 


342 ±8 


73 ±8 


18 ±6 




PN009 






31 ±5 


100 ±5 












PN010 


100 ± 20 




220 ± 10 


650 ± 10 




224 ± 5 


27 ±6 






Planetary nebulae (halo population) 


PNhOOl 






414 ±6 


1240 ± 10 


9±3 


294 ±3 


25 ± 3 






PNh002 




18 ±5 


421 ± 7 


1240 ± 10 




259 ±4 


22 ±5 






PNh003 






470 ± 10 


1380 ± 10 




310 ±10 


40 ± 10 






PNh004 






380 ± 20 


1090 ± 20 




260 ± 20 








PNh005 






120 ± 10 


360 ± 20 




100 ± 30 








PNh006 


90 ± 20 




360 ±5 


1062 ± 9 




291 ±6 


52 ±6 




14 ±4 


PNh007 


107 ±6 




120 ±2 


345 ±3 


59 ±4 


264 ± 5 


164 ±5 






PNh008 


48 ±6 


5±1 


337 ±3 


989 ±9 


12.2 ± 0.9 


284 ± 2 


41 ±1 






PNh009 






130 ± 10 


380 ± 10 




100 ± 20 








PNhOlO 






140 ± 20 


420 ± 10 




100 ± 20 








Unclassified 


X001 












100 ±6 


28 ±6 


20 ±6 




X002 


60 ± 10 










100 ± 10 




22 ± 7 




X003 












100 ± 20 








X004 












100 ± 10 








X005 












100 ±7 


25 ±8 






X006 












100 ± 10 


70 ± 10 






X007 












100 ±4 


30 ±3 


16 ±4 


13 ±4 


X008 


60 ± 10 










100 ±6 


41 ±5 






X009 












100 ± 10 








X010 












100 ± 10 









Note. — This table is organized into sections similarly to Table [T1 Extinction correction has not been applied. Line fluxes arc reported relative to 
H/3 — 100; however, line flux ratios for lines with large wavelength separations arc unreliable for spectra where the ADC was not functioning (sec Table |TT. 
For objects where H/3 is not detected, fluxes are instead normalized relative to another line whose flux is given as 100. Table [2] is published in its entirety 
in the electronic edition of the Astrophysical Journal and at https://www.cfa.harvard.edu/-nsandGrs/papers/M31/summary.html A portion is shown here 
for guidance regarding its form and content. 



20 Sanders et al. 



TABLE 3 

Derived quantities for M31 HII regions 



ID 


A v 


Direct 


Z94 


KD02 


N06 N2 


N06 03N2 


PT05 


PVT ONS 




(mag) 


(log(0/H)+12) 


(log(0/H)+12) 


(log(0/H)+12) 


(log(0/H)+12) 


(log(0/H)+12) 


(log(0/H)+12) 


(log(N/H)+12) 



HII001 






HII002 






HII003 






HII004 






HII005 






HII006 





6 ±0.8 


HII007 






HII008 






HII009 


2 


6 ±0.4 


HII010 





2 ±0.3 


HII011 


1 


1±0.4 


HII012 





6 ±0.5 


HII013 





49 ±0.07 • 


HII014 





7 ±0.8 


HII015 





6 ±0.4 


HII016 






HII017 





5 ±0.6 


HII018 






HII019 






HII020 






HII021 


1 


± 1 


HII022 






HII023 






HII024 






HII025 


2 


1±0.6 


HII026 


1 


4 ±0.8 


HII027 





7 ±0.4 


HII028 


1 


1±0.5 


HII029 





97 ±0.09 


HII030 


1 


1±0.5 


HII031 


2 


± 1 


HII032 





4 ±0.4 


HII033 





6 ±0.5 


HII034 


2 


8 ±0.8 


HII035 





1±0.2 


HII036 


2 


2 ±0.7 


HII037 


1 


0±0.7 


HII038 


1 


4 ±1.0 


HII039 





7 ±0.5 


HII040 


1 


87 ±0.05 


HII041 


3 


1±0.7 


HII042 


2 


0±0.4 


HII043 


1 


5 ±0.9 


HII044 






HII045 





59 ±0.05 8 


HII046 





8 ±0.2 


HII047 








HII048 





33 ±0.04 ■ 


HII049 


1 


2 ±1.0 


HII050 


3 


0±0.6 



33 ±0.04 



9.06 ±0.03 

8.90 ±0.08 
8.87 ±0.09 





8.85 ± 0.05 


8.83 ± 0.06 










8.76 ±0.08 


8.89 ±0.09 










8.68 ±0.09 


8.70 ±0.07 


8.86 ±0.05 






8.72 ± 0.02 


8.74 ± 0.01 


8.63 ± 0.01 


8.63 


8.23 ± 0.02 


7.36 ± 0.02 




8.8 ± 0.1 


9.0 ± 0.2 










8.87 ± 0.07 


9.01 ± 0.09 












9.0 ± 0.1 












8.35 ±0.01 










8.9 ±0.1 


8.9 ±0.2 












9.03 ± 0.01 












9.05 ± 0.03 










8.8 ± 0.1 


8.99 ± 0.10 










8.76 ± 0.08 


8.85 ± 0.09 










8.81 ±0.09 


8.79 ±0.09 


8.93 ±0.07 






8.94 ±0.02 


8.85 ±0.02 


8.78 ±0.02 


8.83 ±0.01 


8.29 ±0.04 


7.46 ±0.03 






9.1 ± 0.1 


8.77 ±0.05 








8.82 ± 0.06 












8.82 ±0.07 


9.1 ± 0.1 










8.6 ±0.2 


8.91 ±0.07 












8.73 ±0.09 










8.9 ±0.1 


9.1 ± 0.1 










8.9 ±0.1 


9.2 ± 0.2 












8.82 ±0.09 


8.85 ±0.07 






8.86 


8.84 


8.85 ±0.01 


8.76 


8.30 ±0.01 


7.47 ±0.01 




8.8 ±0.1 


8.99 ±0.07 












8.8 ±0.1 


8.88 ±0.05 








8.9 ±0.1 










8.68 ±0.01 


8.77 


8.40 


8.39 


8.46 


7.39 ±0.02 


8.96 ±0.04 


8.88 ±0.03 


8.86 ±0.02 


8.86 ±0.01 


8.32 ±0.06 


7.53 ±0.05 




9.18±0.01 










8.96 


8.87 


8.55 


8.62 


8.56 


7.57 ±0.01 




8.7 ±0.2 


8.8 ±0.2 












9.10 ±0.07 









Note. — The abundance diagnostics applied in this table are described in £|2.7I The reported uncertainties are derived by propagation of the line flux uncer- 
tainties, as described in i|2.7| and does not include other systematic effects. When this "statisti cal" uncertainty is less t han 0.01 d ex, we do not repo rt it. Table 
[3] is published in its entirety in the electronic edition of the Astrophysical Journal and at https://www.cfa.harvard.edu/~risanders/papers/M31/summary.html 
A portion is shown here for guidance regarding its form and content. 



Metallicity profile of M31 



TABLE 4 
Derived quantities for M31 

PLANETARY NEBULAE 



ID 






PN001 






PN002 


1 


3 ±0.1 


PN003 





3 ±0.2 


PN004 





09 ± 0.05 


PN005 





4 ±0.5 


PN006 


1 


5 ± 0.3 


PN007 





4 ±0.2 


PN008 





6 ± 0.3 


PN009 






PN010 






PN011 






PN012 





7 ±0.5 




U 


U 


PN014 





01 ±0.09 


PN015 


1 


8 ±0.3 


PN016 





5 ± 0.7 


PN017 





4 ±0.4 


PN018 





0±0.4 


PN019 





6 ± 0.2 


PN020 


2 


0±0.1 


PN021 





1 ± 0.1 


PN022 


1 


0±0.4 


PN023 


2 


± 1 


PN024 


1 


1 ± 0.9 


PN025 






PNhOOl 





1 ± 0.1 


PNh002 






PNh003 





2 ± 0.4 


PNh004 








PNh005 






PNh006 





1 ± 0.2 


PNh007 






PNh008 








PNh009 






PNhOlO 






PNhOll 








PNh012 





17 ±0.07 


PNh013 








PNh014 






PNh015 





3 ± 0.4 


PNh016 





3 ±0.3 


PNh017 





01 ±0.03 


PNh018 





2 ± 0.1 


PNh019 





02 ±0.06 


PNh020 





22 ±0.10 


PNh021 






PNh022 








PNh023 






PNh024 






PNh025 





1 ± 0.2 



log(0/H)+12 



51 ±0.04 



±0.2 



±0.1 



8.61 ±0.07 

8.62 ±0.05 



Note. — 



See caption for Table [3] 



Sanders et al. 



TABLE 5 

Fitted parameters for the abundance profiles of M31 HII regions and PNe 



N Bootstrap Spearman 

Central Slope p p 

(dex) (dex kpc~ 1 ) 



HII regions — All 



A v 


199 


2.31 ±0.26 


-0.0690 ±0.0168 


-0.31 


8 x 10~ 


06 


log([N II] /Ha) 


223 


-0.42 ±0.06 


-0.0087 ± 0.0048 


-0.17 


0.01 




R23 


61 


0.43 ± 0.06 


0.0169 ± 0.0044 


0.54 


8 x KT 


00 


P 


61 


0.26 ± 0.08 


0.0077 ±0.0061 


0.18 


0.18 




log(0/H)+12 (Z94) 


60 


9.10 ±0.06 


-0.0208 ± 0.0048 


-0.57 


2 x 10" 


06 


log(0/H)+12 (KD02) 


136 


8.96 ± 0.06 


-0.0096 ± 0.0049 


-0.33 


1 x 10" 


04 


log(0/H)+12 (N06 N2) 


192 


9.13 ±0.07 


-0.0195 ± 0.0055 


-0.26 


3 x 10- 


04 


log 0/H)+12 (N06 03N2) 


100 


8.98 ±0.08 


-0.0130 ± 0.0068 


-0.19 


0.06 




log(0/H)+12 (PT05) 


18 


8.42 ± 0.09 


-0.0054 ± 0.0064 


-0.08 


0.60 




log(N/H)+12 (PVT ONS) 


52 


7.83 ±0.07 


-0.0303 ± 0.0049 


-0.53 


5 x KT 


05 



HII regions — Stellar 





92 


2.20 ±0.36 


-0.0555 ± 0.0231 


-0.27 


0.01 


log([N II] /Ho) 


98 


-0.34 ± 0.07 


-0.0080 ±0.0056 


-0.14 


0.18 


#23 


14 


0.54 ±0.12 


0.0107 ±0.0096 


0.52 


0.06 


P 


14 


0.04 ±0.11 


0.0218 ±0.0068 


0.55 


0.04 


log(0/H)+12 (Z94) 


14 


9.00 ±0.16 


-0.0152 ± 0.0133 


-0.52 


0.06 


log(0/H)+12 (KD02) 


47 


8.92 ±0.12 


-0.0037 ± 0.0092 


-0.16 


0.27 


log(0/H)+12 (N06 N2) 


72 


9.24 ±0.11 


-0.0224 ± 0.0082 


-0.39 


6 x 10~ 04 


log(0/H)+12 (N06 03N2) 


25 


9.10 ±0.13 


-0.0191 ± 0.0110 


-0.33 


0.10 


log(0/H)+12 (PT05) 


10 


8.18 ±0.19 


0.0077 ± 0.0141 


0.33 


0.35 


log(N/H)+12 (PVT ONS) 


10 


7.80 ±0.16 


-0.0283 ±0.0112 


-0.42 


0.23 



HII regions — Diffuse 



A v 


107 


2.75 ±0.37 


-0.1066 ± 0.0261 


-0.40 


2 x lO - 


05 


log([N II] /Ha) 


125 


-0.37 ±0.06 


-0.0177 ±0.0049 


-0.26 


4 x 10~ 


03 


R23 


17 


0.38 ±0.07 


0.0210 ± 0.0044 


0.52 


2 x 10- 


04 


P 


17 


0.33 ±0.11 


0.0033 ± 0.0089 


0.08 


0.60 




log(0/H)+12 (Z94) 


16 


9.16 ±0.07 


-0.0248 ± 0.0048 


-0.56 


6 x 10- 


or, 


log(0/H)+12 (KD02) 


89 


9.01 ±0.04 


-0.0160 ± 0.0031 


-0.43 


2 x 10- 


(Vy 


log(0/H)+12 (N06 N2) 


120 


9.09 ±0.09 


-0.0195 ± 0.0070 


-0.20 


0.03 




log(0/H)+12 (N06 03N2) 


7.5 


8.93 ±0.11 


-0.0102 ± 0.0095 


-0.11 


0.35 




log(0/H)+12 (PT05) 


38 


8.51 ±0.09 


-0.0113 ±0.0065 


-0.16 


0.33 




log(N/H)+12 (PVT ONS) 


12 


7.83 ±0.08 


-0.0300 ± 0.0058 


-0.54 


2 x 10- 


04 


PNe — Disk 


A v 


333 


0.58 ±0.08 


0.0070 ± 0.0071 


0.04 


0.51 




log([N II]/Hq) 


277 


-0.57 ±0.05 


-0.0023 ± 0.0044 


-0.05 


0.38 




R23 


118 


1.20 ±0.03 


-0.0016 ± 0.0027 


-0.07 


0.42 




P 


118 


0.93 ±0.02 


-0.0027 ±0.0020 


-0.17 


0.04 




log(0/H)+12 (direct) 


51 


8.47 ±0.09 


-0.0056 ± 0.0076 


-0.11 


0.45 





Note. — The bootstrap fitting reported in this table is described in ^|3.2| The central value 
and slope define the best-fit line to the data with respect to the de-projected radius in M31. 



